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The massless Dirac electron is expected to display half-integer quantum Hall effect (QHE) at the filling factor 
u — ±1/2, ±3/2, ±5/2, • • • , reflecting the parity anomaly. However, it has so far been unable to materialize 
such a series due to the fermion doubling problem inherent to the chiral symmetric lattice system. For instance, 
the QHE in graphene displays a series z/ = ±2,±6,±10, - - because of the 4-fold degeneracy of each Landau 
level. We demonstrate that the half-integer series v — ±1/2, ±3/2, ±5/2, • • • can arise in silicene, which is 
a silicon cousin of graphene, when we break the chiral, time-reversal and inversion symmetries all together by 
applying electric field and photo-irradiation simultaneously. A prominent hallmark is the emergence of a single 
Dirac-cone (SDC) state, which contains one massless Dirac cone and three massive Dirac cones in the Brillouin 
zone. A SDC state is a topologically protected semimetal since it emerges along the boundary between two 
distinctive topological insulators. The half-integer QHE is predicted to appear in such a SDC state. 



The study of graphene has begun with the experimental dis- 
covery of an unconventional quantum Hall effect (QHE) with 
Hall plateaux at the filling factor v = ±2, ±6, ±10, • • • dH, 
implying the 4-fold degeneracy of each Landau level. The 
Hall conductivity increases by e'^/h when the Fermi energy 
crosses one Landau level. If there were no degeneracy the 
Hall conductivity would be half-integer quantized jsl], 



ctr = ±(n ± l/2)e'^/h, n = 0, 1, 2, 



(1) 



It reflects the parity anomaly of the massless Dirac fermion JH 
This is also a result of the particle-hole symmetry to- 
gether with the fact that the Landau level of the massless Dirac 
fermion resides at exactly zero energy [|61. However the "half 
integer" is hidden in graphene under the 4-fold degeneracy as- 
sociated with the spin and valley degrees of freedom. Though 
there is a long history in quest for the genuine half-integer 
QHE given by ©, it has never been discovered. This is be- 
cause the necessary condition is the emergence of a single 
massless Dirac-cone state, which we call a single Dirac-cone 
(SDC) state for brevity. The main obstacle is the Nielsen- 
Ninomiya theorem, which requires an even number of mass- 
less Dirac fermions in the chiral symmetric lattice systemd?!]. 
In the surface of the topological insulator, a SDC state is real- 
ized because two Dirac cones are split into the top and bottom 
surfaces However, this resolution is only superficial: 
In a slab geometry electrons on the top and bottom surfaces 
contribute to the Hall conductance all together, yielding the 
integer QHE01 

Silicene iTlU 131 [TtI 155] may well open a way to material- 
ize the genuine half-integer series (1). Silicene consists of a 
honeycomb lattice of silicon atoms with buckled sublattices 
made of A sites and B sites. The states near the Eermi en- 
ergy are tt orbitals residing near the K and K' points at oppo- 
site corners of the hexagonal Brillouin zone. The low-energy 
dynamics in the K and K' valleys is described by the Dirac 
theory as in graphene, but electrons are massive due to a rel- 
atively large spin-orbit (SO) coupling. Silicene is a quantum 
spin Hall (QSH) insulator[14], which is a particular type of a 
two-dimensional topological insulator. 

The Kane-Mele model tl4 ll5] describes well silicene. In 
addition, a staggered potential and the Haldane term^l^ may 








— j^-^ 





FIG. 1 : Illustration of a SDC state, which consists one massless Dirac 
cone and three masive Dirac cones in the valence band. In this in- 
stance, up-spin (red) electrons are massless at the K point but massive 
at the K' point, while down- spin (blue) electrons are massive both at 
the K and K' points. 



be introduced to the system by applying external electric 
field|17] and a circularly polarized photo-irradiation [18, 19], 
respectively. A rich topological phase diagram is obtained 
since the electric field breaks the space inversion symmetry 
and the photo-irradiation breaks the time-reversal symmetry. 
Furthermore, the mass term breaks the chiral symmetry, inval- 
idating the Nielsen-Ninomiya theorem in the system. We are 
thus able to create a SDC state| 19], which contains one mass- 
less Dirac cone and three massive Dirac cones [Fig{T]. In this 
paper we analyze the QHE. We demonstrate that, by tuning 
the electric field and the photo-irradiation, the massless elec- 
tron leads to the half-integer QHE given by Q. The other 
three masses are different, and the plateaux series is predicted 
to be±l/2,±,3/2,±5/2,---. 

We analyze the Kane-Mele-Haldane Hamiltonian in the 
presence of the staggered sublattice potential term. The tight- 
binding model is| 15, 17, 19], 

+ ^7^ Yl ^^i4aCj/3±^Xl^^^^^i«Cm, (2) 

where cj^ creates an electron with spin polarization a at site i, 
and (i, j) / {{ijj)) run over all the nearest/next-nearest neigh- 
bor hopping sites. We explain each term, (i) The first term 
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{a1 ) Band gap ( Aq=0) K 1 K' T ^^^^ ^^"^^ ( An= Aso/2) K ' 



FIG. 2: Phase diagram in the {Ez, An) plane. A circle shows a point 
where the Dirac-cone structure is given in FigQ Heavy lines repre- 
sent phase boundaries indexed by Kry and Sz =ti- There appear a 
SDC state along the line, which is characterized by a single gapless 
Dirac cone at the Kry point with spin Sz. The topological charges 
(C, Cs) are indicated. We give the band gaps along the holizontal 
lines at Aq = and = Aso/2 (dotted line) in Fig|3] 



represents the usual nearest-neighbor hopping with the trans- 
fer energy t = 1.6eV. (ii) The second termlisi |2Q|i repre- 
sents the effective SO coupling with Aso = 3.9meV, where 
a = {(Jx^cFy^ cFz) is the Pauli matrix of spin, with Vij = +1 
if the next-nearest-neighboring hopping is anticlockwise and 
Uij = —1 if it is clockwise with respect to the positive z axis, 
(iii) The third term represents the Haldane interaction induced 
by the photo-irradiation, where Aq = fiVp^^l]"^ with Vt the 
frequency and A the dimensionless intensity 1 18, 19]. (iv) The 
fourth term jivl] is the staggered sublattice potential term. Due 
to the buckled structure the two sublattice planes are separated 
by a distance, which we denote by 21 with I = 0.23 A. It gen- 
erates a staggered sublattice potential oc 2£Ez between silicon 
atoms at A sites and B sites in electric field Ez . 

The low-energy effective Dirac Hamiltonian in the momen- 



tum space reads 11 15 



(3) 

where = = b.h x lO^m/s is the Fermi velocity with 

the lattice constant a = 3. 86 A, and r = (rx^Ty^Tz) is the 
Pauli matrix of the sublattice pseudospin. The factor r] = ±1 
distinguishes the K and K' points. 

The Hamiltonian ([3]) describes a four-component Dirac 
fermion indexed by the spin = ±1 and the valley r] = ±1. 
We use Sz =ti and r] =K,K' for indices to make their mean- 
ing clear. The coefficient of is summarized as A^^ . It is the 
mass of the Dirac electron with the spin Sz in the valley r], 



Al =r]SzXso-^Ez^v^n. 



(4) 



which may be positive, negative or zero. The band gap is 
given by 2 1 A^^ | . It is a remarkable feature that we can control 
the spin- valley dependent mass A^^ by changing the electric 
field Ez and the photo-irradiation strength Aq. 

Silicene is an ideal playground to explore topological in- 
sulators and topological phase transitions[8i,t9J. The topolog- 
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FIG. 3: (al,bl) The band gaps as a function of Ez with the fixed 
values of Aq = and Aq = ^Aso. (a2,a3,b2,b3) The structure of 
four Dirac cones at special points as indicated by circles in Fig|2l 
(a2) The band gap is open in the QSH insulator. (a3) There are two 
massless Dirac cones in the SVPM state. (b2,b3) These are typical 
SDC states. We have set Aso = 0.2t for illustration. 



ical quantum numbers are the Chern number C and the Z2 
index. When the spin Sz is a good quantum number, as is 
the case of the present model, the Z2 index is identical to 
the spin-Chern number Cs modulo 2. They are defined by 
C = + Cf + Cl + Cf and C. = + - Cf - Cf ), 
where C^^ is the summation of the Berry curvature in the mo- 
mentum space over all occupied states of electrons with spin 
5 2 in the Dirac valley It is straightforward to calculate A^^ , 
which are given by llM 11 IH 



C^,^ = ^sgn(A^^ 



(5) 



as a function of the Dirac mass A^ . 

It is a trivial task to derive the topological phase diagram in 
the {Ezj Xq) plane, just by calculating the topological num- 
bers (C, Cs) with the help of ([5]), as illustrated in FigjJl A 
topological phase transition occurs when the sign of the mass 
changes. There are four phase boundaries [FigjS], which are 
given by solving A^^ = 0, corresponding to Sz =ti and 
T] =K,K'. Let us review two typical cases. 

We increase Ez from Ez = while keeping Xq = in 
the phase diagram [Fig 121. We show the band gap 2| A^^ | as a 
function of Ez in FigOal). The Dirac cones at =0 are 
illustrated in FigI3a2). As Ez increases, two Dirac masses 
decrease and vanishes (A^ = A^ = 0) at Ez = Ecr as in 
FigHal). Silicene becomes a spin- valley polarized semimetal 
(SVPM), where the Dirac cones are illustrated in Figj3la3). 
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Doubling of massless fermions occurs. As Ez increases fur- 
ther, the gap opens again and siHcene becomes a trivial bund 
insulator (BI) [FigH. 

We next increase while keeping Xq = Aso/2 in the 
phase diagram [FiglS. Two Dirac masses (A^ and ) de- 
crease and the other two (A^ and A^ ) increase, as illustrated 
inlSbl). We find A^ = at = ^E^,, where up-spin 
electrons are massless at the K point while all others are mas- 
sive, as the Dirac cones show in[3b2). This is a SDC state 
(SDCl). As Ez increases further, the gap opens again and 
silicene becomes a spin-polarized quantum anomalous Hall 
(SQAH) insulator. We also find Af = Sit Ez = f £^cr, 
where down-spin electrons are massless at the K' point while 
all others are massive, as the Dirac cones show in Figj3lb3). 
This is another SDC state (SDC2). 

We proceed to apply a homogeneous magnetic field B = 
V X A = (0, 0, — 5) with B > along the z axis to silicene. 
The Hamiltonian is given by making the minimal substitution, 
hki Pi = hki + eAi, in the Hamiltonian ([3]). We introduce 
a pair of Landau-level ladder operators. 



(a1 ) Landau levels (Aq=0) at K point (b1 ) Landau levels (Aq=Aso/2) at K point 



V2h 



lB{Px-iPv) 
V2h 



(6) 



satisfying [a,a^ = 1, where £3 = \/h/eB is the magnetic 
length. The Hamiltonian iJ^ is block diagonal and given by 



H„ 





HI 



with the diagonal elements being 



Hi = 
Hf = 



(7) 

(8) 
(9) 



in the basis {ipA, "^b} ■ 

It is straightforward to solve the eigen equation of H^^ . The 
lowest Landau level is special, where the eigenvalue is 



El=riAl {Ez). 



(10) 



The eigenstate describes electrons when E'^^ > and holes 
when E'^ < 0. For higher Landau levels they are given by 



±EN = ±^{nwc)^N^{Al)^, Ar = l,2,---, (11) 

where ± for the electron and hole states. 

The present result is applicable also to the QHE in graphene 
by setting the Dirac mass zero for all spin- valley components 
(A^^ = 0). We have the 4-fold degeneracy in all Landau 
levels. Provided the electron-hole symmetry is required, the 
series reads v = ±2,±6,±10,---,asis consistent with the 
experimental factJH 13] • 

We next consider silicene without electric field (Ez = E^x) 
and photo-irradiation (Aq = 0), where the Dirac mass reads 




(a2) Landau levels (An=0) at K' point (b2) Landau levels (Aq=A3o/2) at K' point 




(a3) Hall plateau (£z=£cr,AQ=0) (b3) Hall plateau (£z=£cr,AQ=Aso/2) 
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FIG. 4: (al,a2,bl,b2) Landau levels as a function of Ez at special 
points of Aq = 0, |Aso. The two lines crossing the ^^-axis rep- 
resent the lowest Landau levels. (al,a2) The lowest Landau level 
resides at exactly zero energy for Ez — =b£^cr when \n = 0. 
All levels are doubly degenerate. (a3) Hall plateaux develop at 
h> = zbl,±3,=b5,---. (bl,b2) The lowest Landau level resides at 
exactly zero energy for Ez — =b|^cr, ±f^cr when Aq = |Aso. 
All levels are nondegenerate. (b3) Hall plateaux develop at = 
±1/2, ±3/3, ±5/5, • • • . We have set Aso = 0.2t for illustration. 



A^^ = r]SzXsO' Each levels are doubly degenerate, and hence 
the series reads||22|] = ±1, ±3, ±5, • • • . 

Finally we analyze the SDC state at Ez = ^E^r and 
= ^80. where only electrons with Sz =t at the K point 
are massless. All levels are nondegenerate, and the series read 
ly = ±l/2,±3/3,±5/5,---,by requiring the electron-hole 
symmetry. There are many other SDC states along the phase 
boundary [Fig IS, where the half-integer QHE ([T} appears. 

In conclusion, we have investigated the Kane-Mele- 
Haldane Hamiltonian in the presence of the staggered sublat- 
tice potential term. SDC states emerge along the phase bound- 
aries in the phase diagram [Fig 121 . The SDC state contains 
one massless Dirac cone and three massive Dirac cones. Such 
a state is allowed because we break the chiral, time-reversal 
and inversion symmetries explicitly. We have demonstrated 
the occurrence of half-integer QHE in the SDC state. 

The model ([3]) describes well silicene under electric field 
and photo-irradiation. There are actually small correction 
terms due to Rashba interactions [19J. However we can nu- 
merically show [22] that these terms produce quantitatively 
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minor corrections and not change our conclusions. Further- 
more we can make the following arguments on this point. 

In general, when the topological numbers change across 
the boundary, the band must close and yield gapless modes 
in the interface. Otherwise the topological numbers cannot 
change the values across the interface. Now, a SDC state ap- 
pears along the boundary between two distinctive topological 
insulators. Since they are topologically protected, the SDC 
state is also topologically protected against small perturba- 
tions. Consequently, we may say that the SDC state is a topo- 
logical semimetal. 

One might question the validity of the electron-hole sym- 
metry when higher order corrections are included into the 
model. Even so, we wish to emphasize that the QHE is a 
phenomenon taking place in the vicinity of half- filling, where 
the Dirac system (O describes the system very well and the 
electron-hole symmetry is a good symmetry. Indeed, the ex- 
perimetal observation of the unconventional QHE at u = 
±2,±6,±10,-- - in graphene[l, 2] supports this conclusion 
because higher order corrections are common to graphene and 
silicene. 
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